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ABSTRACT 

We perform ray-tracing simulations evaluating the effect of a cD galaxy on the strong lensing 
properties of five galaxy cluster halos obtained from N-body simulations. The cD galaxy is 
modelled using both axially symmetric and elliptical models and assuming several masses for 
its dark matter halo. The effect of the cD orientation with respect to the mass distribution of 
the host galaxy cluster is also investigated. The simulations are carried out in an open and a 
flat model universe with cosmological constant. We find that the enhancement of the cluster 
lensing cross sections for long and thin arcs due to the presence of a massive cD at the cluster 
centre is typically less than 100%, depending on the model used for the cD galaxy and its 
orientation. The impact of the cD on the cluster efficiency for producing radially magnified 
images is larger only for those clusters whose lensing cross section for radial arcs is very small 
in absence of the central galaxy. We conclude that the presence of a cD galaxy at the cluster 
centre can only moderately influence the cluster efficiency for strong lensing and in particular 
fails to explain the discrepancy between the observed number of giant arcs in galaxy clusters 
and their abundance predicted from lensing simulations in the currently most favoured ACDM 
model. 



1 INTRODUCTION 



Many authors have pointed out that the observed abundance of 
gravitational arcs whose length-to-width ratio exceeds a given 
threshold can be used for constraining the cosmological model and 
its parameters. Indeed, the occurrence of those rare events caused 
by highly non-linear effects in cluster cores is determined by sev- 
eral factors related to the geometry of the universe, as well as to the 
abundance, evolution and internal structure of the lensing clusters, 
which all depend on cosmology. 

Investigating the lensing properties of a large sample of nu- 
merically modelled galaxy clusters using the ray-tracing tech- 
nique, Bartelmann et al. (1998) showed that the expected num- 
ber of giant arcs, usually defined as arcs with length-to-width 
ratio exceeding 10 and apparent B-magnitude less than 22.5 
(Wu & Hammer 1993 1, changes by orders of magnitude between 
different cosmological models. In particular, they found that the 
expected number of giant arcs is smaller by a factor of ten in a flat, 
low density universe with a matter density parameter £Iq = 0.3 and 
a cosmological constant = 0.7 (ACDM) than in an open, low- 
density cosmological model with £^o = 0-3 and = (OCDM). 

Observations of the abundance of gravitational arcs in galaxy 
clusters seem to be consistent only with the predictions for 
an open universe (Bartelmann et al. 1998 Le Fevre et al. 19941 
Gioia & Luppino 1994| |Luppino et al. 1999| iGladders et al. 20031 
Zaritsky & Gonzalez 2003 1. This is in pronounced disagree- 
ment with other observational results, in particular those 
obtained from the recent experiments on the cosmic mi- 
crowave background l de Bernardis et al. 2002 1 and the observa- 
tions of high-redshift type-la supernovae l Perlmutter et al. 19971 



IPerlmutter 19 981 which all suggest instead that the cosmological 
model most favoured by the data is spatially flat and dominated by 
a cosmological constant. 

Several studies tried to resolve this inconsistency, 
but satisfactory explanations remain to be found. Simi- 
lar estimates of the number of long and thin arcs were 
attempted adopting isothermal analytic instead of numeri- 
cal cluster lens models jCooray, Quashnock & Miller 1999| 
IKaufmarm"& Straumann 2000}. The results of these studies differ 
from those of Bartelmann et al. (1998) in that they are almost 
insensitive to the cosmological constant and predict similar 
abundances of arcs in the ACDM and in the OCDM models. 
However, Meneghetti, Bartelmann & Moscardini (2002) recently 
showed that analytic models are inadequate for quantitative studies 
of arc statistics, because asymmetries and substructures in the 
deflector play a relevant role for determining the galaxy cluster 
efficiency for strong lensing and cannot properly be taken into 
account when using simple analytic models for modelling cluster 
lenses. Moreover, Bartelmann et al. (2003) demonstrated that 
elliptical cluster models with the density profile found in numerical 
simulations by Navarro et al. (1996) were able to reproduce the 
relative change of an order of magnitude in the arc cross sections 
produced by clusters in the ACDM and OCDM models. 

By investigating the properties of giant gravitational arcs in 
the Red-Sequence Cluster Survey, Gladders et al. (2003) found a 
surprisingly high frequency of arc systems in high redshift clusters 
(z > 0.64). Their determination of lensing incidence is in agreement 
with the results by Zaritsky & Gonzalez (2003), who measured the 
abundance of strong gravitational lensing events in the Las Cam- 
panas Distant Cluster Survey. As the authors of these studies sug- 
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gest, some property associated with clusters at early times might re- 
sult in a higher efficiency for strong lensing. Torri et al. (2003) have 
recently studied the effects on strong cluster-lensing cross sections 
induced by major merger events within the clusters. Their results 
indicate that the cross section for giant arcs can change by almost 
one order of magnitude while a massive substructure crosses the 
cluster centre. Numerical simulations show that mergers are fre- 
quent events in clusters. Therefore, they may offer an explanation 
as to why so many arc systems are observed in both high and low 
redshift clusters. Further investigations into the impact of mergers 
on arc statistics are now in progress. 

Meneghetti et al. (2001) also investigated the effect of clus- 
ter galaxies on arc statistics but found that they do not introduce 
perturbations strong enough for significantly changing the number 
of arcs and the distributions of lengths, widths, curvature radii and 
length- to- width ratios of long arcs. 

However, the possible effects due to the presence of a massive 
cD galaxy near the cluster centre were not investigated there. The 
centres of massive galaxy clusters are generally dominated by such 
very massive (~ 1O 13 M0) galaxies, which could in principle no- 
ticeably affect the strong-lensing properties of their host clusters. 
In fact, due to their more concentrated dark matter halos, they may 
steepen the inner slope of the cluster density profile and push the 
cluster critical curves to larger distances from the centre. Thus, the 
length of the critical curves may be increased, and thus the proba- 
bility for long arcs to form. Moreover, cD galaxies may help their 
host clusters to reach the critical central surface density for produc- 
ing critical curves and becoming efficient strong lenses. 

Understanding the effect of massive central galaxies on the 
strong lensing properties of galaxy clusters is an important issue 
also because several studies attempted using arc statistics for con- 
straining the density profiles of lenses or the nature of the dark mat- 
ter composing their halos <Wu & Mao 1996 Molikawa et al. 19991 
Oguri et al. 200T1 IMolikawa & Hattori 20011 |Oguri et al. 2002) 
Meneghetti et al. 20019. The conclusions from all these studies 
may change should the effect of a central massive galaxy turn out 
to be relevant. 

Thus, we investigate in this paper whether the presence of a 
massive cD galaxy at the centre of a galaxy cluster can signifi- 
cantly change the cluster's ability to produce strong lensing events 
and reconcile the lensing cross section of ACDM clusters with the 
observations. For doing so, we study the lensing properties of a 
sample of five galaxy clusters numerically simulated in both the 
ACDM and the OCDM cosmological models, and we measure their 
efficiency for producing tangential and radial arcs before and after 
the inclusion of a cD galaxy. The central galaxy is modelled using 
both axially symmetric and elliptical models and assuming a range 
of virial masses and possible orientations with respect to the mass 
distribution of the host cluster. 

The plan of the paper is as follows. In Section|2| we use sim- 
ple analytic models for qualitatively estimating the impact of a cD 
galaxy on the strong lensing properties of galaxy clusters. In Sec- 
tion |3| we describe the numerical models used and the methods 
adopted for the lensing simulations. The results are discussed in 
Section|4| Finally, conclusions are summarised in Section|5| 



2 EXPECTATIONS 

We begin by investigating with the help of analytical models what 
impact we can expect from cD galaxies on the strong-lensing prop- 
erties of galaxy clusters. We will use two different analytic lens 



models, one with a singular isothermal density profile, and the other 
with the density profile found in numerical simulations by Navarro 
et al. (1996; hereafter NFW). 



2.1 Axially symmetric models 

The radial density profile of a singular isothermal sphere (hereafter 
SIS) is given by 
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where a v is the velocity dispersion of the halo. The NFW density 
profile is instead given by 



p(r) = 
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where p s and r s are characteristic density and distance scales, re- 
spectively (see Navarro et al. 1997). These two parameters are not 
independent. The ratio between r s and the radius r2oo within which 
the mean halo density is 200 times the critical density is called 
concentration, c = r2oo />s. As results from numerical simulations 
show, the concentration parameter c can be expressed as a func- 
tion of the halo virial mass, which thus is the only free parameter. 
The concentration also depends on the cosmological model, im- 
plying that the lensing properties of haloes with identical mass are 
different in different cosmological models if they are modelled as 
NFW spheres. The NFW halo profile falls off steeper than isother- 
mal at radii beyond r s , but flattens towards the halo centre. These 
different features lead to markedly different lensing properties of 
the NFW compared to the SIS model (see the discussions in Per- 
rotta et al. 2002 and Meneghetti et al. 2003). 

Several algorithms were developed for computing the concen- 
tration parameter from the halo virial mass <Navarro et al. 19971 
IBullock et al. 2001llEke et al. 20011 . They are all based on the find- 
ings that numerically simulated haloes tend to be more concen- 
trated the earlier they form, and that their central density reflects the 
mean cosmic density at their formation time. Since haloes of lower 
mass form earlier in hierarchical models than haloes of higher 
mass, the concentration is a decreasing function of the halo mass. 

In this work, we adopt the algorithm proposed by Navarro et 
al. (1997), which first assigns to a halo of mass M a collapse red- 
shift z co n defined as the redshift at which half of the final mass is 
contained in progenitors more massive than a fraction / of the final 
mass. Then the density scale of the halo is assumed to be some fac- 
tor C times the mean cosmic density at the collapse redshift. They 
recommend setting / = 0.01 and C = 3 x 10 3 because their numer- 
ically determined halo concentrations were well fit adopting these 
values. 

In the case of axially symmetric models the computation of the 
deflection angles reduces to a one-dimensional problem. We define 
the optical axis as the straight line passing through the observer and 
the lens centre and introduce the physical distances perpendicular 
to the optical axis on the lens and source planes, t, and T|, respec- 
tively. We then fix a length scale £o on the lens plane and define the 
dimensionless distance x = from the lens centre. By project- 
ing £o on the source plane, we define a corresponding length scale 
T|o = (D s /D\) on the source plane, where D s and D\ are the an- 
gular diameter distances to the source and lens planes, respectively. 
Like on the lens plane, we define a dimensionless distance from the 
optical axis y = r\/r\o on the source plane. 

Using this dimensionless formalism, the lensing potential of a 
SIS lens can be written as 
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if cjg = 47t(o",,/c) 2 {D\D\ S /Ds) is chosen as a length scale, where 
D\ s is the angular-diameter distance between lens and source. The 
derivative of the lensing potential with respect to x is the reduced 
deflection angle at distance x, 
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For an NFW sphere, we choose cjo = r s an d define K s = 
Ps^Zct 1 , where Z cr = [c 2 /(AnG)\ [£> s /(AAs)] is the critical sur- 
face mass density for strong lensing. Its lensing potential then reads 



V(x)=4K s g(x) 
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jMeneghetti et al. 2003|IBartelmann 199S1 . 
2.2 Elliptical model 

Lensing by galaxy clusters can only crudely be described by axially 
symmetric models because their generally high level of asymmetry 
and substructure changes their lensing properties qualitatively and 
substantially. We thus perturb the axially symmetric lens models 
such that their lensing potentials and ^5) acquire elliptical iso- 
contours. We define the ellipticity as e = 1—b/a, where a and b are 
the major and minor axes of the ellipse, respectively, and introduce 
it into J5} by substituting 
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where x\ and x 2 are the two Cartesian components of x, x 2 = x\- 
x\. The Cartesian components of the new deflection angle are 
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where &(x) is the unperturbed (i.e. axially-symmetric) deflection 
angle at the distance X from the lens centre. 

2.3 Estimates 

Arcs form close to the tangential critical curve of a lens from 
sources close to the tangential caustic. The size of the lensing cross 
sections for the formation of tangential arcs is strictly connected to 
the extent of the tangential critical curves and caustics: the longer 
the critical curves and the corresponding caustics are, the larger are 
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Figure 1. The solid and the dashed curves show the ratios between the tan- 
gential critical radii of a cluster of virial mass 7.5 x 10 14 /r'M with and 
without a cD galaxy of virial mass 10 13 h~ [ Mq added to its centre. They 
were obtained modelling the cluster with a SIS or NFW density profile, 
respectively. The dotted curve shows the ratios between the radial critical 
radii, obtained modelling the cluster with an NFW density profile. The ra- 
tios are given as functions of the ellipticity e of the cluster potential, as 
described in the text. Since critical curves are not circular for e > 0, we use 
their maximum cluster-centric distance instead. 



the lensing cross sections. For axially symmetric lens models, the 
tangential critical curve encloses a circle which encloses an average 
convergence of unity, so that 
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defines the radius t of the tangential critical curve. For singular 
isothermal spheres, 

A ' 
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which is approximately proportional to the lens mass. A similar 
closed expression cannot be given for an NFW sphere. The tan- 
gential critical radius of a set of concentrically superposed SIS 
lens models is simply the sum of the tangential critical radii of 
the individual components. This implies that adding a cD galaxy 
of ~ 1O 13 /i _1 M extends the tangential critical curve of a cluster 
of ~ 5 x 10 14 h~ l Mq by a few per cent only. 

However, the situation changes remarkably if the cD galaxy, 
modelled as a SIS, is added to a cluster with an NFW density pro- 
file. The comparatively flatter density profile near the cluster centre 
causes the cD galaxy to have a much stronger effect, as Fig.Qillus- 
trates. 

The solid and dashed curves shown in the figure are the ratios 
between the tangential critical radii of a cluster with virial mass 
7.5 x 10 A _1 Af0 with and without a cD galaxy of virial mass 
10 13 h~ l Mq added to its centre. They are obtained for clusters with 
SIS and NFW density profiles, respectively. The curves show how 
these ratios change with the ellipticity e of the cluster's lensing po- 
tential. The tangential critical curves are circles only for e = 0; for 
e > 0, we use their maximum cluster-centric distance instead. 

Two results can be read off Fig. Q First, the effect of a cD 
galaxy of fixed mass on the tangential critical curves of a cluster 
of fixed mass is much more pronounced if the cluster has an NFW 
rather than a SIS density profile; and second, the effect of the cD 
decreases as the ellipticity of the cluster increases because the in- 
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Figure 2. The curves show, as functions of lens redshift, the maximum 
extent of tangential critical curves and corresponding caustics of a clus- 
ter with NFW density profile and virial mass 7.5 x lO 14 /r'A/ , with and 
without a cD galaxy with singular isothermal density profile and virial mass 
10 13 hT l Mq. An ellipticity of e = 0.3 was assumed for the cluster's lens- 
ing potential, and the source redshift is Zs — 1 • Obviously, the impact of a 
cD galaxy is highly insensitive to the lens redshift. The curves end where 
the maximum extent of the tangential critical curves drops below 10" be- 
cause we are interested in the formation of large arcs only, which require a 
minimum size of the tangential critical curves. 

creased shear of the cluster makes the additional surface mass den- 
sity relatively less important. For the axially symmetric NFW clus- 
ter, the cD extends the tangential critical curve by as much as 38%, 
although it has only a few per cent of the cluster's mass! Numerical 
simulations suggest ellipticities of ~ 0.3, for which the effect drops 
~ 12%. 

Since we have good reasons to believe that clusters are better 
described by NFW than by SIS density profiles, the figure illus- 
trates that an estimate of the effect of cD galaxies on strong lensing 
by clusters based on SIS models alone can be highly misleading. 
In addition, the change with cluster ellipticity is important because 
it shows that the enhanced gravitational shear due to asymmetries 
and substructure in the cluster potential reduces the impact of the 
cD galaxy. 

It is a second interesting question whether the impact of a 
cD galaxy on a cluster lens changes with cluster redshift. Again, 
for two concentrically superposed SIS lenses, the relative change 
in the tangential critical radii is exactly independent of redshift, 
but this does not need to hold true for cluster lenses with NFW 
density profile. We show in Fig. |2| as a function of lens redshift 
the maximum extent of the tangential critical curve and the corre- 
sponding caustic for a cluster with NFW density profile and virial 
mass 7.5 x W li h~ l M@, with and without a cD galaxy with sin- 
gular isothermal density profile and virial mass \Q n hT x M & . The 
cluster's lensing potential is assumed to have ellipticity e = 0.3, in 
good agreement with values found in numerical simulations. 

The curves end where the maximum extent of the tangential 
critical curve drops below 10", because arcs formed near smaller 
tangential critical curves can hardly be classified as large arcs. The 
curves illustrate that the relative change in the dimensions of tan- 
gential critical curves and caustics is remarkably insensitive to the 
lens redshift, even if the cluster is assumed to have an NFW density 
profile. 

Finally, we consider the effect of a cD galaxy on a cluster's ef- 
ficiency for producing radial arcs. The dotted curve shown in Fig.Q 



is the ratio between the radial critical radii of a cluster with and 
without a cD galaxy added to its center. The virial masses of the 
cluster and the cD galaxy are the same as in the previously dis- 
cussed tests. Again, the cD galaxy is modelled as a SIS while the 
cluster has an NFW density profile. As the figure shows, the impact 
of the cD galaxy is negligible when the cluster is axially symmet- 
ric. Indeed, the radial critical curve of such a cluster is a very small 
circle around the cluster center. Adding a cD galaxy modelled as a 
SIS does not change much because it dominates the density profile 
only close to the cluster core. Since for a SIS the derivative of the 
deflection angle, doc/dx, is zero everywhere, while the condition 
for having a radial critical curve in an axially symmetric model is 
da/dx = 1, the extent of the radial critical curve is not changed 
by the central galaxy. For it to be efficient, the radial critical curve 
or parts of it need to be pushed outward where the density and 
deflection-angle profiles are no longer dominated by the cD galaxy, 
and once this is achieved by the ellipticity, the cD galaxy can in- 
troduce considerable change. We thus expect that the impact of cD 
galaxies on the efficiency for producing radial arcs of realistic clus- 
ter models might be much larger than for tangential arcs. 

We now turn to numerical cluster models in order to test the 
impact of cD galaxies on cluster lensing with more realistic cluster 
mass distributions. 



3 NUMERICAL MODELS 
3.1 Simulated Clusters 

We investigate the lensing properties of a sample of five numeri- 
cally simulated cluster-sized dark matter halos kindly made avail- 
able by the GIF collaboration l Kauffmann et al. 1999 1. The same 
clusters were used by Bartelmann et al. (1998). They were obtained 
from A^-body simulations performed in the framework of several 
cosmological models. In this paper, we only use the simulations 
for a flat model with cosmological constant (ACDM), which has 
£Iq = 0.3, = 0.7 and h = 0.7, and in an open model without 
cosmological constant (OCDM, Cl = 0.3, h = 0.7). 

The initial matter density in these models is per- 
turbed about the mean according to a CDM power spectrum 
IBond & Efstathiou 1984 1 with primordial spectral index n = 1, 
normalised such that the local abundance of massive galaxy 
clusters is reproduced (e.g. Viana & Liddle 1996). The complete 
list of cosmological parameters in these simulations is given in 
Bartelmann et al. (1998). 

For each cosmological model, an initial simulation with A' = 
256 3 particles in a box of 141 h _i Mpc side length was run. The 
mass of individual particles is 1.4 x 10 10 h~* M@. Clusters were 
obtained from initial cosmological simulations as follows. High 
density regions were identified by means of a standard friends-of- 
friends algorithm, selecting only the dense cores of the collapsed 
objects. Around the centres of those, all particles were collected 
which lie within a sphere of radius 1.5/j~' Mpc, corresponding to 
the Abell radius. These objects are taken as clusters. For our anal- 
ysis, we took the five most massive clusters in the simulation vol- 
umes. 

As shown by Bartelmann et al. (1998), the maximum effi- 
ciency for the production of long and thin arcs is reached by these 
clusters when their redshift is in the range 0.2<z<0.4. Since our 
analytic estimates showed that the lens redshift is highly irrelevant 
for the impact of a cD on strong cluster lensing, we only take the 
simulation snapshots at redshift zl = 0.275, where the number of 
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simulated arcs is largest and thus the uncertainties in the numeri- 
cally determined cross sections are smallest. 

3.2 Lensing simulations 

The lensing properties of the numerical clusters are studied using 
the ray-tracing technique. Our method has been fully described in 
several earlier publications (see Meneghetti et al. 2000, 2001). We 
thus give only a brief description here, referring the reader to those 
papers for further detail. 

Each numerical cluster is centred within a cube of 3h~ Mpc 
side length. Within this box, we compute the three-dimensional 
density field p by interpolating the mass density on a regular grid 
of 256 3 cells, using the Triangular Shape Cloud method (TSC; 
see Hockney & Eastwood 1988). Finally, we project p along the 
three independent box sides, obtaining three surface density maps 
for each cluster. These are used as lens planes. Adopting this thin 
screen approximation is justified because the distances between the 
sources and the lens and between the lens and the observer are 
much larger than the physical sizes of the galaxy clusters. 

We then propagate a bundle of 2048 x 2048 light rays through 
the central quarter of each of these maps. We focus on the cen- 
tral parts of the clusters because our analysis concerns their strong- 
lensing properties only. Thus, our ray-tracing simulations resolve 
scales of order ~ 0.2", at the lens redshift. The deflection angles of 
each ray are computed by directly summing the contributions from 
each mass element of the cluster. 

Once the deflection angles are known, we reconstruct the im- 
ages of a large number of extended background sources. They are 
all placed on the same plane, located at the redshift z s = 1. Al- 
though real sources are distributed in redshift, putting them at a 
single redshift is justified for the present purposes because the lens 
convergence changes very little with source redshift if the lens red- 
shift is substantially smaller, as is the case here. 

The sources have elliptical shape, with axial ratios randomly 
drawn from the range [0.5, 1], and their area is equal to that of a 
circle of 1" diameter. We initially distribute them on a coarse reg- 
ular grid in the source plane. Then, their spatial density is itera- 
tively increased towards the caustic curves, where the lensing ef- 
fects strengthen. This increases the probability of producing long 
arcs and thus the numerical efficiency of the method. In order to 
compensate for this artificial source-density increase, we assign in 
the following statistical analysis to each image a statistical weight 
proportional to the area of the grid cell on which the source was 
placed. 

Image reconstruction and classification is done following the 
technique introduced by Bartelmann & Weiss (1994) and used in 
Meneghetti et al. (2000, 2001, 2003). For each image, we measure 
its length and width. This results in a catalogue of simulated images 
which is subsequently analysed statistically. 

3.3 Inclusion of the cD galaxy 

Being linear function of mass, the total deflection angle of a ray 
passing through a mass distribution is the sum of the contributions 
from each mass element of the deflector. Therefore, in the case of a 
galaxy cluster, we can decompose the cluster lens into its smoothed 
dark matter component, plus the granular component contributed 
by its galaxy population (see also Meneghetti et al. 2000). For both 
the cluster and the galaxies, the main constituent is given by the 
dark matter which forms their halos. Our model of the cluster con- 
taining a cD galaxy can thus be fairly simple; we take the smoothed 
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dark matter distribution obtained from the numerical simulations 
described above, and introduce a dark-matter halo resembling the 
galaxy. For each ray traced through the lens plane, we compute the 
deflection angle by summing the contributions from the cluster it- 
self and the galaxy haloes. 

Using this approach, the reduced deflection angle of each ray 
parameterised by its cell numbers (i,j) on the grid in the image 
plane is given by 

a, 7 = a^. + a^ D , (13) 

where a?j and a?? are contributions to the deflection angle from 
the original cluster and the cD galaxy, respectively. 

The values of af? depend on the model for the cD galaxy. We 
first apply two axially symmetric models, namely spheres with the 
NFW or the singular isothermal density profile. Second, we also ap- 
ply the pseudo-elliptical NFW lens model l Meneghetti et al. 2003) 
in order to account for the possible elongation of the matter distri- 
bution of the cD galaxy. 

cD galaxies typically appear to be of elliptical shape, with 
isophotal axis ratios b/a ~ 0.8 ( Porte r"et al . 1 99 it . Moreover, the 
orientation of the brightest cluster ellipticals is usually not random, 
but correlates well with that of their host cluster iSastry 1968 
|Rood & Sastry 1972| lAustin & Peach 19741 Icarter 1980 
Bing geli 1982| IStruble & Peebles 19851 |Rhee & Katgert 1987 
Lambas et al. 1988 Garijo et al. 1997). Asymmetries in the lens- 
ing matter distribution are known to improve the ability of the 
cluster to produce long and thin arcs. We also expect that the 
impact of a cD galaxy described by an elliptical model is largest 
when its orientation is aligned with the elongation direction of the 
host galaxy cluster. In the case of different orientations, the cD 
galaxy tends to circularise the mass distribution of the cluster in its 
central region. In order to quantify this effect, we carry out two sets 
of simulations; in the first, we randomly choose the orientation of 
the cD galaxy inside the cluster, while in the second the orientation 
is chosen such that the directions of the major and minor axes 
of the galaxy align with the major and minor eigenvalues of the 
cluster's deflection angle field, respectively. The galaxy ellipticity 
is assumed to be e — 0.2 l Porter et al. 199TY 

The position of the cD galaxy is chosen to coincide with the 
minimum of the deflection angle field of the numerical galaxy clus- 
ter. This choice intentionally maximises the contribution of the cD 
galaxy to the total deflection angles. 



4 RESULTS 

4.1 Convergence maps and critical curves 

We describe in this Section how the cD galaxy changes the pro- 
jected mass distribution of the numerical clusters. The surface den- 
sity map of a lens can be reconstructed from the deflection angles 
once the effect of the cD galaxy has been included. Indeed, the 
scaled surface density K (or convergence) at each position x on the 
lens plane is given by 

K(x) = 5M = ^V-a(f). (14) 

Some examples for convergence maps are shown in Fig. 
The original convergence map of the cluster, i.e. before adding the 
cD galaxy, is plotted in the top left panel. Then, starting from the 
top central panel, we show the convergence maps for all cD galaxy 
models we use, followed by the convergence map of the cluster in- 
cluding the cD galaxy: the NFW sphere (top central and top right 
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Figure 3. Examples of convergence maps for one of the clusters in the GIF sample and for the cD galaxy models. Starting from the top-left panel: (a) cluster 
without the cD galaxy; (b) cD galaxy of mass M c d = 5 x 'A) ' li '.I/, with NFW density profile; (c) cluster with the cD galaxy modelled as in (b); (d) cD 
galaxy with the same mass, but with SIS density profile; (e) cluster with the cD galaxy modelled as in (d); (f) pseudo-elliptical NFW model randomly oriented 
with respect to the cluster mass distribution; (g) cluster with the cD galaxy modelled as in (f); (h) pseudo-elliptical NFW model with its orientation aligned 
with that of the cluster mass distribution; (i) cluster with the cD galaxy modelled as in (h). The scale of each panel is ~ 200/i~'kpc or ~ 70 arcsec. 



panels), SIS (middle left and middle central panels), the pseudo- 
elliptical NFW model both randomly oriented (middle right and 
bottom left panels) and aligned with the orientation of the clus- 
ter mass distribution (bottom central and bottom left panels). The 
virial mass of the cD galaxy is the largest we have tested in our 
simulations (5 x 10 13 bT x Mr.)). The angular scale of each plot is 
approximately ~ 50" . 

The contours of the convergence maps show that the presence 
of the cD galaxy affects the matter distribution of the cluster only in 
its very central part. The size of the region where the convergence 
enhancement is significant depends on the model which is used to 
describe the cD galaxy. For example, the SIS has a steeper density 
profile at the centre and thus its contribution to the host cluster's 



surface density is significant only in a very limited region. On the 
other hand, due to its shallower density profile, the mass of a cD 
galaxy modelled as an NFW sphere or as a pseudo-elliptical NFW 
model is distributed across a larger region around the cluster centre. 
Moreover, when the pseudo-elliptical model is used, convergence 
contours appear to be more extended and particularly elongated if 
the cD galaxy is aligned with the orientation of the underlying clus- 
ter mass distribution. 

Arcs form along the critical curves, where the determinant 
of the Jacobian matrix A(x) = dy/dx vanishes. Given that the 
dimension-less coordinates on the source plane y and on the lens 
plane 1 are related by the lens equation, y = x — 6t(x), the position 
of the critical curves can be easily found through numerical deriva- 
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Figure 4. Examples of critical curves and caustics. The first row of panels shows how the critical curves for one projection of the most massive cluster 
(M ~ 10 h~ M@) in the GIF sample change in response to the inclusion of the cD galaxy. In each panel, we draw the original critical curves (bold curves) 
and those modified by the cD (thin curves). We show the results obtained by modelling the cD galaxy as an NFW sphere (first column from the left); using 
the pseudo-elliptical model randomly oriented (second column) and aligned with the cluster mass distribution (third column), and as a SIS (fourth column). 
The panels in the second row show the corresponding caustic curves. Similarly, we plot in the third and in the fourth row of panels the critical curves and the 
caustics for one projection of the least massive cluster in the sample (M ~ 3 x 10 14 h~ [ Mq). The mass of the cD galaxy is M c d = 5 x 10 i3 /i~'Mq, and the 
background cosmology is the ACDM model. The physical scale of each panel is ~ 265 h~ ' kpc. 



tives of the deflection angles. Then, using the lens equation, we can 
find the position of the caustics, near which sources must be located 
in order to be strongly magnified. 

We show in Fig. some examples of critical curves and caus- 
tics for the most massive (first two rows of panels) and the least 
massive (last two rows of panels) clusters in our sample (in the 
ACDM model). In each panel, we plot the results before adding 
the cD galaxy (thick curves), and after including the contribution 
to the deflection angles from a cD galaxy with a virial mass of 
5 x 10 I3 h~ Mr., (thin curves). In the first column on the left, we 
show the results obtained by modelling the cD as an NFW sphere; 
in the second column, we plot the critical curves and caustics found 
modelling the galaxy as pseudo-elliptical, randomly oriented NFW 
models; in the third column, we show the results for the pseudo- 
elliptical NFW model aligned with the orientation of the underly- 
ing cluster matter distribution; and in the fourth column we finally 
show the critical curves and the caustics obtained modelling the cD 
as a SIS. 



As expected, the critical curves and the caustics appear to be 
pushed towards larger cluster-centric radii when the cD galaxy is 
included into the simulations. The largest effect is produced by the 
pseudo-elliptical cD if its orientation is aligned with that of the host 
cluster. The singular isothermal cD has the least effect on the crit- 
ical curves. Of course, the impact of the cD galaxy is different on 
clusters with different masses. For example, the caustics of the most 
massive lens demonstrate that they are affected by the presence of 
the cD only far from the cusps. This is not the case for the least mas- 
sive cluster, where also the cusp positions are pushed away from the 
lens centre. Also, the critical curves appear to be more extended 
once the cD is added to the least massive cluster. Of course, this is 
due to the fact that, by adding the cD, we are changing the cluster 
mass in its central region by a different percentage. Similar results 
are found for clusters in the OCDM model. 
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Figure 5. Relative change in the cross sections for arcs with length-to-width ratios exceeding 5 (left column), 7.5 (central column) and 10 (right column) 
as a function of cluster mass for the numerically simulated clusters in the ACDM model. Results are shown for three different masses of the cD galaxy: 
5 x 10 /; \l . (top panels), I0 i! /; M. (middle panels) and 5 x 10 13 h M@ (bottom panels). Filled circles, triangles and squares mark the results obtained 
modelling the cD as an NFW sphere, a pseudo-elliptical NFW model with random orientation, and aligned with the orientation of the host cluster, respectively; 
open circles show the results found modelling the cD galaxy as a SIS. 



4.2 Tangential arcs 

Tangential arcs are selected from the catalogues of simulated im- 
ages for each numerical cluster model by requiring a minimal 
length-to-width ratio. We then quantify the cluster's efficiency for 
producing this type of arcs by measuring its lensing cross section. 
This procedure is carried out for all numerical clusters before and 
after the inclusion of the cD galaxy. 

By definition, the lensing cross section is the area on the 
source plane where a source must be located in order to be imaged 
as an arc with the specified property. As explained in Sect. 13,21 each 
source is taken to represent a fraction of the source plane. We assign 
a statistical weight of unity to the sources which are placed on the 
sub-grid with the highest resolution. These cells have area A. The 
lensing cross section is then measured by counting the statistical 
weights of the sources whose images satisfy a specified property. If 
a source has multiple images with the required characteristics, its 
statistical weight is multiplied by the number of such images. Thus, 
the formula for computing cross sections for arcs with a property p 
is 



O p =A^w,«; , (15) 

i 

where n; is the number of images of the i-th source satisfying the 
required conditions, and w, is the statistical weight of the source. 

We compute the lensing cross sections for arcs whose length- 
to-width ratio (L/W) exceeds a lower threshold (L/W)m„. The re- 
sults obtained for each cluster with and without the cD galaxy are 
compared by computing the relative change of the cross sections 
due to the presence of the cD. 

We show in Figs. J5J and JSJ the relative change of the lens- 
ing cross sections for long and thin arcs as a function of the cluster 
virial mass for the simulations in the ACDM and in the OCDM 
models. The panels in the left, central and right columns show 
the relative change of the cross sections for arcs with length-to- 
width ratios (L/W) > 5, (L/W) > 7.5 and (L/W) > 10, respec- 
tively. In each panel, we plot the results for all four models used 
to describe the cD galaxy. The top, middle and bottom panels re- 
fer to simulations in which the virial mass of the cD galaxy is 
(5 x 10 12 , 10°, 5 x M.. respectively. 

As expected, the largest variations of the cross sections are 
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Figure 6. Relative change in the cross sections for arcs with length-to-width ratios exceeding 5 (left column), 7.5 (central column) and 10 (right column) 
as a function of cluster mass for the numerically simulated clusters in the OCDM model. Results are shown for three different masses of the cD galaxy: 
5 x 1O 12 /z~'M0 (top panels), 10 h~ Mq (middle panels) and 5 x 10 3 h .W . (bottom panels). Circles, squares and triangles identify different cD galaxy 
models with the same symbols as in Fig. |3J- 



typically found if the cD galaxy is modelled as a pseudo-elliptical 
NFW model whose orientation is aligned with that the host cluster. 
On the other hand, cD galaxies with SIS profiles change the ability 
of the numerical clusters for producing long and thin arcs only by 
a very small amount. 

In the ACDM model, a cD galaxy with mass between M c d = 
5 x 10 12 bT 1 Mq : , and M cD = 10 13 h" 1 M produces maximal varia- 
tions on the order of 40% — 50%. More massive cDs have a stronger 
impact on the strong-lensing efficiency of the clusters: a galaxy 
with mass M C £> = 5 x 10 13 h~ M@ changes the lensing cross sec- 
tion by a maximum amount between 60% and 200%, depending on 
the total cluster mass. The impact of the cD is generally larger in 
the less massive clusters. 

A similar trend is found in the OCDM model, but the vari- 
ations of the cross sections are smaller in this case. For example, 
in the simulations with the most massive cDs, the cross sections 
change by approximately 40% — 80% only. This behaviour was 
expected because the clusters in the OCDM model are generally 
more compact compared than those in the ACDM model. Includ- 
ing the cD, the mass in the very central part of the clusters changes 
less compared to the clusters in the ACDM model. Cross sections 



for arcs with different minimal length-to- width ratios show similar 
variations. 

4.3 Radial arcs 

The appearance and location of radial arcs within the lensing clus- 
ters is determined by several factors, among them the slope of the 
projected density profile and the central density of the lens. The 
steeper the density profile is, the closer to the centre the radial arcs 
move. Moreover, the higher the central surface density is, the larger 
is the extent of the radial critical line. 

Radial arcs are quite rare events. Many clusters in our sam- 
ple exhibit very short radial critical curves and are thus not very 
efficient in producing radially distorted images. In particular, the 
numerical clusters in our sample have significantly smaller cross 
sections for radial arcs in the ACDM than in the OCDM model. 
The reason is that the clusters in the ACDM model are less concen- 
trated than in the OCDM model. Given that we increase the mass 
in the inner part of the numerical clusters by quite a large amount 
by adding a massive cD galaxy, we expect that the capability of the 
cluster models for producing radial arcs should be increased when 
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Figure 7. Relative change of the cross section for radial arcs versus original cross sections (i.e. before the inclusion of a cD galaxy). The left and right panels 
show the results obtained in the ACDM and the OCDM models, respectively. We use circles, squares and triangles to identify different cD galaxy models, 
using the same symbols as in Figs. J5J and 



a cD galaxy is included. For investigating this question, we again 
compare the lensing cross sections for radial arcs of the numerical 
models before and after including a cD galaxy. 

Radial arcs are identified from the complete sample of dis- 
torted images using the technique described by Meneghetti et 
al. (2001). It consists in selecting those arcs for which the measured 
radial magnification at their position exceeds a lower threshold. The 
relative change of the lensing cross section for this type of arcs as 
a function of the original lensing cross section of clusters not con- 
taining a cD galaxy is shown in Fig. J7}. The cD galaxy mass is 
5 x 10 13 h~ l Mq. In the left and right panels, we plot the results for 
the cluster models in the ACDM and in the OCDM models, respec- 
tively. Again we use circles, squares and triangles for identifying 
different cD galaxy models, using the same symbols as in Figs. J5J 
and 0. 

As the figure shows, the cD galaxy increases the capability 
of the cluster models for producing radial arcs, but this enhance- 
ment depends on the original size of their lensing cross sections. 
Clusters which were less efficient in producing radial arcs increase 
their lensing cross sections by roughly 500%, while this increment 
is smaller for lenses which already produced a higher number of 
radial arcs. In particular, the effect of the cD galaxy on the lensing 
cross sections of clusters is significantly larger in the ACDM model 
compared to the OCDM model. One of the clusters in our sample 
(the least massive one) did not produce any radial arcs before a 
cD galaxy was included, thus the results for only four clusters are 
shown in Fig. Q. 



5 SUMMARY AND CONCLUSIONS 

We investigated the effect of massive cD galaxies on the strong- 
lensing properties of galaxy clusters. The main motivation behind 
our study is to identify possible reasons for the pronounced fail- 
ure of numerical cluster-lensing models in the strongly favoured 
ACDM cosmology in reproducing the observed abundance of giant 
luminous arcs. 

Although inadequate for quantitatively reliable results, our an- 



alytic model demonstrated three interesting aspects of the problem. 
First, the effect of cD galaxies is much stronger in clusters with 
NFW rather than singular isothermal density profile because of the 
flatter central density profile of the former. Second, cD galaxies 
are expected to be much less efficient in asymmetric than in axi- 
ally symmetric clusters because the stronger gravitational tidal field 
in asymmetric clusters makes the relative contribution of the addi- 
tional central mass component less important. Third, the impact of 
cD galaxies on the strong-lensing properties of clusters is expected 
to be almost independent of lens redshift. 

We then carried out ray-tracing simulations using massive nu- 
merically simulated galaxy clusters in the ACDM and OCDM cos- 
mologies to which we added cD galaxies in four different ways; 
either modelled as singular isothermal spheres or modelled with 
the NFW density profile, either axially symmetric or elliptically 
distorted; and if distorted, either oriented randomly or aligned with 
the orientation of the cluster mass distribution. All cD galaxies are 
placed at the minima of the deflection-angle maps of the cluster 
models. They have masses of 5 x 10 12 , 10 13 , or 5 x 10 3 ft M@, 
which are added to the cluster mass. 

Note that these choices yield conservative results because they 
tend to exaggerate the possible effects. The cD-galaxy masses are 
very high, in particular compared to the least massive cluster of our 
sample, and they are added to the clusters without correspondingly 
decreasing the cluster mass. 

Although we find relative enhancements of large-arc cross 
sections by up to ~ 200% for our least massive cluster if a cD 
galaxy with 5 x 10 13 hr 1 M and an elliptically distorted NFW den- 
sity profile is added with its orientation aligned with the mass dis- 
tribution of the cluster, such a situation can hardly be considered 
realistic because that cD galaxy has ~ 20% of the cluster mass. In 
more realistic cases of cD galaxies with < 10 13 h^ 1 Mq, tangential- 
arc cross sections are increased by not more than ~ 50%. It should 
be noted, however, that we have to compare the relative increase 
between clusters in ACDM and OCDM models, which is typically 
not more than ~ 30%. 

Cross sections for the formation of radial arcs are much more 
affected because they are highly sensitive to the exact central den- 
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sity profile of the cluster lenses. The presence of a cD galaxy can 
multiply the radial-arc cross sections by factors of a few, typically 
up to five in the ACDM model. This reinforces the finding by other 
authors that radial-arc statistics potentially is a highly significant 
tracer for central density profiles of clusters, but as such less suit- 
able for constraints on cosmological parameters. Currently, also, 
the data base for radial arc statistics is rather too poor for any mean- 
ingful conclusions. In fact only few radial arcs have been discov- 
ered so far and all of them have been found in clusters which are 
dominated by a cD galaxy (e.g. AC 114, MS 2137; Natarajan et 
al. 1998, Hammer et al. 1997) or by one or two giant ellipticals 
(e.g. A 370; Bezecourt et al. 1999). 

We thus conclude from our conservative estimates of the im- 
pact of cD galaxies on strong-lensing cross sections by galaxy 
clusters that they may increase the arc-formation probability by 
perhaps up to ~ 50% in realistic situations, but certainly by far 
not enough for explaining the discrepancy between simulations in 
ACDM models and the observed abundance of arcs. 

This result supports the view that we are still a long way from 
understanding what exactly makes real clusters efficient strong 
lenses. Numerous observational effects need to be taken into ac- 
count in addition to careful numerical models before detailed com- 
parisons between observed and simulated arcs can be carried out. 
The recent suggestion by Wambsganss et al. (2003) that arc statis- 
tics in ACDM models can be reconciled with the observed fre- 
quency of arcs by adopting a broader source redshift distribution 
is certainly no solution either because the arc statistics problem as 
posed by Bartelmann et al. (1998) concerns the ability of a well- 
defined sample of massive clusters to form arcs from a photometri- 
cally limited sample of sources, rather than the strong-lensing op- 
tical depth of a section of the Universe for arbitrarily faint sources. 
It is very likely that careful comparisons between numerically sim- 
ulated and observed strong-lensing clusters will tell us much about 
the detailed mass distribution and the history of mass assembly in 
individual clusters. 
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